Qubit decoherence by low-frequency noise 
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We have derived explicit non-perturbative expression for decoherence of quantum oscillations in 
a qubit by low-frequency noise. Decoherence strength is controlled by the noise spectral density at 
zero frequency while the noise correlation time r determines the time t of crossover from the 
to the exponential suppression of coherence. We also performed Monte Carlo simulations of qubit 
dynamics with noise which agree with the analytical results and show that most of the conclusions 
are valid for both Gaussian and non-Gaussian noise. 
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Despite the large number of successful demonstrations 
of coherent quantum oscillations in individual 0, |3j S 
Si Hi IE Qi 111 and coupled 9] Josephson-junction qubits, 
quantitative understanding of these oscillations is so far 
limited. The main area of discrepancy between experi- 
ment and theory is qubit decoherence. The typical qual- 
ity factors of reported oscillations, while not as large as 
required by potential applications in quantum computa- 
tion, are still quite large in physics term (typically not 
less than 20 30) . This fact should imply weak decoher- 
ence describable by the standard perturbation theory in 
qubit-environment coupling (see, e.g., 0|). Several basic 
features of this theory, however, do not agree with exper- 
imental observations. Most importantly, observed decay 
time T2 of coherent oscillations is typically shorter than 
the energy relaxation time Ti even at optimal qubit bias 
points 0, 0, llH where perturbation theory predicts no 
pure dephasing terms. Another discrepancy is between 
the observed two-qubit decoherence rate ^ and its val- 
ues that can be obtained from the perturbation theory 
under natural assumptions 12]. 

Qualitatively, the basic reason for discrepancy between 
Ti and T2 is the low-frequency noise that can reduce 
T2 without changing significantly the relaxation rates. 
Mechanisms of low-frequency, or specifically 1//, noise 
exist in all solid-state qubits: background charge fluc- 
tuations for charge-based qubits impurity spins or 
trapped fluxes for magnetic qubits Manifestations 
of this noise are observed in the echo-type experiments 
Low-frequency noise for qubits is also created by 
the electromagnetic fluctuations in filtered control lines. 

The goal of our work is to develop quantitative theory 
of low-frequency decoherence by studying qubit dynam- 
ics under the influence of noise with small characteristic 
amplitude vq and long correlation time r. In the case 
of Gaussian noise, we obtained explicit non-perturbative 
expression describing decay in time of coherent qubit os- 
cillations. The strength of decoherence in this expression 
is controlled by the noise spectral density at zero fre- 
quency, VqT. For long correlation times r 3> A"-'^, where 
A is the qubit tunnel amplitude, the spectral density VqT 
can be large even for weak noise vo ^ A and our analyt- 



ical results are exact as function of VqT in this limit. We 
also performed direct numerical simulations of the low- 
frequency qubit decoherence. The simulation results con- 
firm analytical expressions for Gaussian noise and show 
that our main conclusions: cross-over from the l/Vt to 
the exponential suppression of coherence at time t ~ t; 
and the strength of decoherence controlled by the noise 
spectral density Sy{0) at zero frequency remain valid for 
non-Gaussian noise. The decoherence rates for Gaussian 
and non-Gaussian noise behave, however, quite differ- 
ently as functions of 5^(0). In physics terms, the two 
noise models we study correspond, respectively, to low- 
frequency electromagnetic noise, and noise of localized 
excitations in the situation when only a few excitations 
with comparable time scales are coupled to the qubit. 

The Hamiltonian of a qubit with a fluctuating bias 
energy v{t) (see inset in Fig. lb) is: 

H = -^[Aa, + ie + v{t))a,], (1) 

where e is the average bias, and u's here and below de- 
note Pauli matrices. In this work, we mostly focus on 
the situation when the noise v{t) has characteristic cor- 
relation time r, i.e., the noise correlation function and its 
spectral density can be taken as 

{vit)v{t')) = i^o^eH*'*'!/^ , S^ic.) = , (2) 

where vo is the typical noise amplitude and (...) denotes 
average over different realizations of noise. We assume 
that the temperature T of the noise-producing environ- 
ment is large on the scale of the cut-off frequency l/r, 
and it can be treated as classical. (In the regime of in- 
terest, l/r <C A, the temperature can obviously be still 
small on the qubit energy scale.) 

The two effects of the weak noise on the dynamics 
of the qubit ^ are the transitions between two energy 
eigenstates with energies ±17/2, = (A^ -I- e^)^/^, and 
"pure" (unrelated to transitions) dephasing that sup- 
presses coherence between these states. Within the stan- 
dard perturbation theory, the transition rate is propor- 
tional to Sv{^) — Vq/^I'^t. One can see that the condi- 
tion of weak noise vq ^ A makes the transition rate small 
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compared both to A and 1/r ensuring that the perturba- 
tion theory is sufficient for the description of transitions. 
As discussed quahtatively in the introduction, the fact 
that the noise correlation time is long, r ^ A~^, makes 
the perturbation theory inadequate for the description of 
pure dephasing. For low-frequency noise, a proper (non- 
perturbative in Vqt) description is obtained by looking 
at the accumulation of the noise-induced phase between 
the two instantaneous energy eigenstates. If vq <C A, one 
can determine the rate of accumulation of this phase by 
expanding the energies in noise amplitude v{t). Also, in 
this case the dephasing rate is larger than the transition 
rate and can be calculated disregarding the transitions. 
The factor F(t) describing suppression in time of coher- 
ence between the two states (i.e., suppression of the off- 
diagonal element pi2 of the qubit density matrix in the 
energy basis: puit) — i^(i)pi2(0)e~'^*) can be written 
then as follows: 



F{t) = (exp{~i 



ev{t') A2v2(<') 



dt'}) . (3) 



For Gaussian noise, the correlation function ^ deter- 
mines the noise statistics completely. In this case, it is 
convenient to take the average in Eq. (PJ by writing it 
as a functional integral over noise. For this purpose, and 
also for use in the numerical simulations, we start with 
the "transition" prohahility p{vi,V2, St) for the noise to 
have the value V2 a time dt after it had the value vi: 



p{vi,V2,t) = [2^«o'(l - e^2**/^)]-i/2 X 
^ 1 iv2 - v,e-''/-)\ 



(4) 



Using this expression we introduce the probability of spe- 
cific noise realization as pQ{vi)-p{vi, V2, Sti)-p{v2, fs, (5^2) • 
where po(w) = {2ttvq)~^^^ expl—v"^ /2vq} is the sta- 
tionary Gaussian probability distribution of v. Taking 
the limit 6tj ^ we see that the average over the noise 
can be written as the following function integral: 



X exp 



v{o)^ + v{ty 

4vi 
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dv{Q)dv{t)Dv{t') 



dt'{T'v'+v') \ . (5) 



Since the average in Eq. (0) with the weight (|3J) is now 
given by the Gaussian integral, it can be calculated 
straightforwardly: 



F(t) =Fo(t)exp 
Fo{t) = e*/^''[coah{vt/T) 



I ft , vt . 

^ ' 2 coth h ly] 

T 2t 



(6) 
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aiBhivt/T)]-^/^ 



where = -^/l + 2ivQA'^T/fl^ and a = ervo/Qi''^^^. 

Equation ^ is our main analytical result for dephasing 
by the Gaussian noise. To analyze its implications, we 
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FIG. 1: The rate 7 of exponential qubit decoherence at long 
times t ^ T for e = and (a) Gaussian and (b) a model of 
the non- Gaussian noise with characteristic amplitude wo and 
correlation time r. Solid lines give analytical results: Eq. © 
in (a) and Eq. I13II in (b). Symbols show 7 extracted from 
Monte Carlo simulations of qubit dynamics. Note different 
scales for 7 in parts (a) and (b). Inset in (b) shows schematic 
diagram of qubit basis states fluctuating under the influence 
of noise v{t). 



start with the case e = 0, where pure qubit dephasing 
vanishes in the standard perturbation theory. Dephasing 
© is still non- vanishing and its strength depends on the 
noise spectral density at zero frequency 5*1,(0) = 2vqT 
through = Vl + is, s = 3^(0)/ A. For small and large 
times t Eq. simplifies to: 



F{t) 



l+t/T 



where 



^=27 



l + t/T + ist/2T 
2V^e-('^+*'5)7(l -1 



(l + s2)l/2 + l 



1/2 



1/2 



1 



(7) 



(8) 



Besides suppressing coherence, the noise also shifts the 
frequency of qubit oscillations. The corresponding fre- 
quency renormalization is well defined for i 3> t: 

1/2 

(9) 
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Suppression of coherence ^ ior t <^ t can be qual- 
itatively understood as the result of averaging over the 
static distribution of noise v. In contrast to this, at large 
times t ^ T, the noise appears to be J-correlated, the 
fact that naturally leads to the exponential decay iQ. 
This interpretation means that the two regimes of decay 
should be generic to different models of the low-frequency 
noise. Indeed, they exist for the non-Gaussian noise con- 
sidered below, and are also found for Gaussian noise with 
1/f spectrum Crossover between the two regimes 

takes place at t ~ r, and the absolute value of F{t) in 
the crossover region can be estimated as (1 + s"^)^^/^, 
i.e. s determines the amount of coherence left to decay 
exponentially. The rate © of exponential decay shows 
a transition from the quadtratic to square-root behav- 
ior as a function of Sv{0) that can be seen in Fig. la, 
which also shows the decay rate extracted from numeri- 
cal simulations of Gaussian noise. (Numerical procedure 
is discussed below.) One can see that our analytical and 
numerical results agree well for quite large noise ampli- 
tude Vq. 

Non-zero qubit bias e leads to additional dephasing 
F{t)/Fo{t) described by the last exponential factor in Eq. 
© . One can see that similar to Fq (t) additional dephas- 
ing exhibits the crossover at t ~ r from "inhomogenious 
broadening" (averaging over the static distribution of the 
noise v) to exponential decay at i ^ r. In contrast to 
FQ{t), the short-time decay is now Gaussian: 



ln[ 



Foity 



v'^Tt/{i + is{A/ny) , t^T. 



We see that, again, the rate of exponential decay de- 
pends non-trivially on the noise spectral density 5*1,(0), 
changing from direct to inverse proportionality to S'„(0) 
at small and large s, respectively. 

Our approach can be used to calculate the rate of ex- 
ponential decay at large times t for Gaussian noise with 
arbitrary spectral density S'^(w). Such a noise can be 
represented as a sum of noises (0 and appropriate trans- 
formation of variables in this sum enables one to write 
the average over the noise as a functional integral similar 
to (|SJ. For calculation of the relaxation rate at large i, 
the boundary terms in the integral |(SJ can be neglected 
and it is dominated by the contribution from the "bulk" 
which can be conveniently written in terms of the Fourier 
components 



Then, (...) = /i?z;...exp{-(l/2)EJi;„|V^.K)}- 
Gombining this equation and Eq. ^ we get at large t: 



F(t) = cxp^-^ 



(iwln(l + iS'^(w)AVf^^) 



(10) 



For unbiased qubit, e = 0, this equation coincides with 
the one obtained by more involved diagrammatic pertur- 
bation theory in quadratic coupling 
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FIG. 2: The profile of coherent quantum oscillations in an 
unbiased qubit dephased by the non-Gaussian noise with 
characteristic amplitude vq = 0.15A and correlation time 
r — 300A~^ obtained by direct simulation of qubit dynamics 
with noise. Solid line is the exponential fit of the oscillation 
amplitude at large times. Dashed line is the initial 1 / \/t decay 
caused by effectively static distribution of v. 

To check how well the analytic theory described above 
works for finite noise amplitude wq, and to see how sen- 
sitive the results are to the assumption of the Gaussian 
noise, we performed Monte Garlo simulations of the qubit 
oscillations under the influence of Gaussian and non- 
Gaussian noise. We looked specifically at the coherent 
oscillations of a qubit with Hamiltonian ^ that start in 
one of the eigenstates of the Cz operator, focusing on the 
case e = 0. The qubit density matrix was averaged over 
up to 10^ realizations of noise. In the case of Gaussian 
noise, realizations were built using the transition prob- 
ability jffl. For non-Gaussian noise we used the model 
of Ref. lla, which should provide an appropriate descrip- 
tion of the situation when a qubit is coupled to several 
fluctuators with similar characteristic time scale r of the 
fluctuations. In this model the fluctuators create random 
qubit bias v which remains constant for some (random) 
time interval after which it is updated and the new value 
remains constant during the next time interval, etc. The 
time intervals between bias updates are taken to be dis- 
tributed according to the Poisson distribution with char- 
acteristic time r. For more direct comparison with the 
Gaussian noise, we assumed Gaussian distribution pq{v) 
of V. The correlation function of v(t) defined in this way 
is given by the same Eq. 

Example of the oscillations dephased by such a noise 
is given in Fig. 2. It shows real part of the off-diagonal 
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element pi2 (t) of the qubit density matrix in the energy 
eigenstates basis. For oscillations starting in one of the 
az eigenstates, pi2{0) = 1/2. Similarly to the case of 
Gaussian noise, we consider only weak noise, Vq ^ A. 
In this case, there is a crossover at t ~ r in the os- 
cillation amplitude from the initial l/\/t suppression of 
coherence due to averaging over static potential distri- 
bution: pi2(t) = pi2{0)/{l + iv^t/Ay/^ (neglecting all 
terms of order t/r), to exponential suppression at < 3> t. 

The rate of the exponential decay can be found analyt- 
ically as follows. Expansion of the average qubit density 
matrix p{t) in the number of "jumps" of v{t) leads to the 
Dyson-like equation for its evolution [l^ : 

p(t)=e-*/^(5(t,0)p(0)5t(i,0)) 
/■* fit' 

+ / ^e-(*-*)/-(5(i,i')p(i')^^(^,0>, (11) 

Jo T 

where (...) denotes the average over the distribution of 
V. For weak noise, introducing slowly-varying amplitude 
r of pi2{t) = r{t)e^^^*, one can reduce Eq. l(TT|l to the 
equation for r(t) neglecting rapidly oscillating terms: 

r(t) =e-*/"r(0)(l+^^'^</A)-l/2 

10 T 



-1/2 



(12) 



With the exponential ansatz for r{t): r{t) cx e '^^ ^)*/'^ 
Eq. H12() gives then equation for the parameter A: 



A = 



dxe 



1 



2A 



^1/2 



(13) 



(Omission of the first term in Eq. (|12|l is justified by 
the final result for the oscillation decay rate 7 = (1 — 
Re[A])/r.) Asymptotics of 7 found from Eq. H13|) are: 



1 

7 — — X 

T 



s78, s<l 
~167r/s2, s> 1 



(14) 



The rate 7(5) evaluated from Eq. H13II is shown in Fig. 
lb, together with the pure dephasing rates found numer- 
ically by fitting the oscillation amplitude (similar to that 
shown in Fig. 2) at t ^ t and subtracting the contri- 
bution 5*^ (A) to dephasing from real transitions. One 
can see that Eq. ifT^ indeed gives an accurate descrip- 
tion of pure dephasing rates. Similarly to the case of 
Gaussian noise, 7 depends only on 5*^(0). In both sit- 
uations, 7 oc Wqt/A^ for small Vqt/A, the fact that 
can be explained by the lowest-order perturbation the- 
ory in qubit energy fluctuations. In the non-perturbative 
regime, however, the behavior of 7 as function of Sv{0) 
is model-dependent and varies from saturation l|14() to 
A/S',„(0)-growth (©. 

In summary, we developed non-perturbative theory 
of qubit dephasing within two models of Gaussian and 
non-Gaussian low-frequency noise and performed Monte 



Garlo simulations of qubit dynamics within these mod- 
els. The theory agrees well with simulations and shows 
that the decoherence strength is controlled by the noise 
spectral density at zero frequency. It allows for general- 
izations in several experimentally-relevant directions and 
should be useful for analysis of observed shapes of quan- 
tum qubit oscillations. 
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